Truncated Levy Flights and Weak Ergodicity Breaking in the 
Hamiltonian Mean Field Model 



&0 
3 



A. Figueiredo, T. M. Rocha Filho, and M. A. Amato 
Instituto de Fisica and International Center for Condensed Matter Physics 
Universidade de Brasilia, CP: 04455, 70919-970 - Brasilia, Brazil 

Z. T. Oliveira Jr 

Departamento de Ciencias Exatas e Tecnologicas, 
Universidade Estadual de Santa Cruz, 
CEP 45662-900 - Brazil and 
International Center for Condensed Matter Physics 
Universidade de Brasilia, CP: 04455, 70919-970 - Brasilia, Brazil 

R. Matsushita 
Departamento de Estatistica 
Universidade de Brasilia, CP: 044 55, 
70919-970 - Brasilia, Brazil 

Abstract 

The dynamics of the Hamiltonian mean field model is studied in the context of continuous time 
random walks. We show that the sojourn times in cells in the momentum space are well described 
by a Levy truncated distribution. Consequently the system in weakly non-ergodic for long times 
that diverge with the number of particles. For a finite number of particles ergodicity is only 
attained for very long times both at thermodynamical equilibrium and at quasi-stationary out of 
equilibrium states. 
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Ergodicity is an important concept in statistical mechanics and essentially states 

that ensemble averages are equal to the time average over a single trajectory of the system or 
that the sojourn time on a given region is proportional to its ensemble measure. For particles 
in an external potential, ergodicity implies that time averages over a single particle trajectory 
is equal to the average over many particles at a fixed time. The latter case has recently been 
observed experimentally for the diffusion of molecules on a nanostructured porous glass j^]. 
Strong ergodicity breaking occurs if some region of the phase space is not accessible by the 
system trajectory. On the other hand weak non-ergodic behavior corresponds to a situation 
where every state can be reached but the occupation statistics is not equal to the ensemble 
measure 4|. The statistical mechanics of weakly non-ergodic systems in the context of 
continuous time random walk (CTRW) was addressed by Rebenshok and Barkai [5]. In this 
approach the system can be in M different states, such that a given observable O admits 
the respective M values O k for k = 1, . . . , M. The time average of this observable is then 

j M M 

O = - — ^tfcOfc, itot = (1) 
hot k= i i=k 

where t k is the residence time, i. e. the total time spent by the system in state k and t tot the 

Is) 

total observation time. The sojourn time t k j is the time spent in state k during the j'-th 

(s) 

visitation, and therefore t k = J2jtkj- Supposing that the sojourn times are independent 
identically distributed (iid) variables and owing to Levy generalized central limit theorem, 
the probability distribution of residence time t k can be described by a one-sided Levy dis- 
tribution fjf^ (t) with Laplace transform 5), 0] : 

fl a \t) e~ st dt = e~^ a/2 , (2) 



for 0<a<2,a = 2 corresponding to the Gaussian distribution, and 7^ is a constant factor. 
The exponent in the right-hand side of eq. ((2]) was chosen to coincide with the exponent 
in the characteristic function of the Levy distribution [6]. An important feature is that 
the moments < > of the Levy distribution diverge for u > a. The distribution of the 
possible different values of the time averages O is given by [5j: 



F^iO) = -- limlmg%^^±^. (3) 



In the limit a -> 2 eq. © reduces to F^(t) = 8(0 - (£>)), where (O) = EklkO k , and 
the coefficients j k are equilibrium (stationary) probabilities. As a consequence ensemble 



and time averages are equal and the system is ergodic. The case of non iid sojourn times 
statistics have been addressed in Ref. |7|. 

In the present work we discuss a system with sojourn time statistics satisfying a truncated 
Levy flight defined as 8|: 



f {a) (t) 



cfW(t), if \t\<L, 
0, otherwise, 



(4) 



where c is a normalization constant and L > 0. Instead of an abrupt truncation at t — L it 
is also possible to use a smooth truncation with a function decaying faster than a power law. 
Such functions have finite moments, and therefore the sum of random variables drawn from it 
converge to a Gaussian distribution. Nevertheless this convergence is extremely slow and for 
a finite sum = Ylk=i ^ w ith drawn from the distribution in eq. (jlj), the distribution 
remains close to a true stable Levy distribution sometimes for large n, displaying a typical 
scaling behavior before a crossover to a Gaussian behavior. The probability density for X^ n ' 
being zero is given for small n by 8]: 

r(i/a) 



P (n) (0) = . . 

with 7 a scale factor and for large n we have a Gaussian distribution: 



(5) 



P (n) (0) = \aV2^n 1/2 



-i 



(6) 



where a is the standard deviation. The crossover value of n from a pure Levy to the Gaussian 
behavior predicted by the central limit theorem increases with L. The statistics of sojourn 
times is a one-sided Levy distribution, but this analysis can be performed by symmetrizing 
the distribution as ff. a \—t) = fj. (t)- The distributions obtained from this sum all collapse 



in the same function under the scaling: 



X -> n- 1/a X and P {n) (X) -> n 1/a P (n) (X). 



(7) 



Examples of weakly non-ergodic systems include laser cooling of tra pped atoms [9(], dif- 

lll . Il2| and glass dy- 



fusion of lipid granules in living cells [10], blinking quantum dots 

namics Q|. In the present letter we study the ergodic properties of a classical Hamiltonian 
system with a long-range interaction considered as a CTRW as a preliminary step of a more 
thorough investigation of more general classes of long-range interacting systems, that have 



drawn some attention in the last two decades 13Ml6l|. and can have some unusual properties 



such as anomalous diffusion, aging, negative heat capacity at equilibrium, non-Gaussian 
quasi-stationary states and a relaxation time to equilibrium diverging with the number N 
of particles. They are characterized by a pair-interaction potential decaying asymptotically 
as r _a , a < D with D the spatial dimensions [l^]. Some important physical systems in 
this category are self-gravitating systems and non-neutral plasmas (see [l4J for more exam- 
ples). Despite their interest their dynamics is not completely understood due to inherent 
difficulties and a few toy models were introduced in the literature in order to simplify the 
understanding of their intricate behavior. Among them we cite the Hamiltonian mean field 



(HMF) model 18j, which has become a sort of ground test for many numerical and analytical 



studies, and defined by the Hamiltonian: 

1 N 1 N 

H(P, 6) = - £P- + — £ [1 - oos(0< - B S )\ . (8) 

i=l M=l 

This system can be interpreted as consisting of N classical rotors globally coupled and with 
unit moment of inertia and its thermodynamical equilibrium properties can be determined 
analytically It has a second order phase transition from a spatially homogeneous to 

an inhomogeneous phase, a rich structure of non-equilibrium phase transitions [20| and the 
great advantage that molecular dynamics simulations scales with the number N of particles 
instead to the usual N 2 scaling allowing simulations with a large number of particles. 
For classical Hamiltonian systems with long-range interactions in the limit iV — > oo 



the one-particle distribution function exactly satisfies the mean field Vlasov equation [22 ]. 
Thence we only have to consider the dynamics of a single particle evolving in the mean field 
of the remaining particles. Quasi-stationary states thus correspond to the infinite number 
of stable stationary states of the Vlasov equation. For finite N small collisional corrections 
must be considered resulting in a secular evolution of the quasi-stationary states. 

The ergodic behavior of the HMF model has been previously discussed by some of the 
authors in Ref. {2l| showing that the time average of the velocity of a single particle is 
equal to the ensemble average over the N particles only after a very long time of the order 
of the relaxation of the system to true thermodynamical equilibrium. Here we discuss this 
behavior in the framework of the CTRW in the momentum space by performing molecular 
dynamics simulations for the Hamiltonian in eq. jSJ). The momentum space is divided in 
cells of width Ap, and each cell is then taken as a different discrete state. We consider 



f(t) 




FIG. 1. Distribution f(t) of sojourn times at the thermodynamic equilibrium for the cells in 
momentum space in the intervals (I) [0.0,0.4] (II) [0.4,0.8], (III) [0.8,1.2] and (IV) [1.2,1.6] for 
N = 100,000 and £ tot = 10 7 . 

here two types of statistically stationary states: thermodynamical equilibrium and a stable 
non-Gaussian (quasi) stationary state with a waterbag one-particle momentum distribution 
p{p) = l/2p if —po < P < Po and zero otherwise. We consider at this point spatially 
homogeneous distributions resulting in a vanishing force in the mean-field limit but with 
small fluctuations due to collisional corrections for finite N which cause small cumulative 
changes in the momenta. The temperature for the equilibrium case is chosen as T = 0.8 
and for the non-equilibrium case we chose po ~ 2.68 which yields an total energy per 
particle e = 0.8. We start by considering the statistics of sojourn times for different cells of 
width Ap = 0.4 shown in Fig. [TJfor the equilibrium case and N = 100,000 particles. The 
distributions obtained are equal up to the noise level, and by supposing that they are also 
independent we have an iid process. The time evolution of the velocity of a single particle 
is shown in the left panel of Fig. [2] for N = 500, 000 particles such that the fluctuations of 
the field are already quite small illustrating clearly the trapping of the velocity in regions 
of the momentum space with rapid movements between different traps. The corresponding 
sojourn times for all cells with width A = 0.4 are shown in the right panel of Fig. |2j The 
right panel of Fig. [2] show the sojourn times in cells of width Ap = 0.4. As the statistics of 
sojourn times are the same the sojourn times are taken from all cells which allows a better 
statistical accuracy in our analysis. 

In order to show that the sojourn time statistics satisfy a truncated Levy distribution 
we perform a sum of n sojourn times randomly chosen among those obtained from our 
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FIG. 2. Left Panel: Velocity of a particle as a function of time for N = 500, 000 and total simulation 
time tf = 10 7 for the Gaussian distribution. Right panel: Corresponding sojourn times. 




FIG. 3. Left Panel: Probability density P$\o) function of the number of summed variables 
for different values of the number N of particles and normalized with respect to the probability at 
zero of the original distribution (n = 1) for the equilibrium state. The straight line is the curve 
obtained from a fitting for the initial part of the curve for N = 200, 00 with a = 0.603. Right 
Panel: The same as the left panel but for the water-bag state, the straight line corresponds to 
a = 0.51. In both panels all curves tend to a t 2 power law with increasing n. 

simulations and also randomly choosing its sign (this is equivalent to symmetrizing the 
distribution). The left panel of Figure [3] shows the probability at zero p(«)(0) as a function 
of n, where the crossover from the scalings in eq. (J5J) to eq. (jH]) is clearly visible. The right 
panel shows the same probability for the non-equilibrium stationary waterbag state. The 
value of the coefficient a in the truncated Levy distribution in ecjs. (12]H| can be determined 



from the sojourn times using maximum likelihood method 



23 



24| and are shown in table |T] 
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10,000 0.6730 0.9993 

50,000 0.5849 0.6868 
100,000 0.5648 0.6451 
200,000 0.5528 0.6129 
500,000 0.5596 0.6265 

TABLE I. Values of a in eq. ([2]) for the thermodynamical equilibrium (a eq ) and waterbag distri- 
bution (a w b). 



for different number of particles. They are obtained by considering the distribution f( a \T) 
as a true (non-truncated) Levy distribution, which is reasonable if L is large enough, but 
can introduce some errors in the estimation of the exponent a. This explains the difference 
between the values of a obtained by this method and those determined using a fitting 
procedure as described in the caption of Fig. [3j 

The second moments < t 2 >= (l/S)Y^k=i^k °f ^ ne distribution of sojourn times, where 
S is the number of sojourn times, are shown for different values of N in figure HI such that 
a power law increase with N is evident, which shows that the effective truncation value 
L diverges with N. Figure [5] shows the tail of the distribution for iV = 50, 000 with the 
straight line corresponding to a power law with a given in table HI Figure E] shows 

the symmetrized distribution functions of sojourn times for a few values of n (left panel) 
and the collapsed distribution obtained from the scaling in eq. (J7|) (right panel), showing a 
good agreement as expected for a Levy distribution. 

From our results we can assert that the HMF model becomes ergodic only after a very 
long time such that the sum of sojourn times obeying a truncated Levy distribution converge 
very slowly to a Gaussian distribution [8] . This sum yields the distribution of residence times 
which at it turn determine the ergodic properties of the system. This time corresponds to the 
crossover between Levy and Gaussian behaviors in the graphics in Figure [31 and for shorter 
times, which are effectively still very long, the system is weakly non-ergodic. Our results 
also show that this crossover time diverges with the number of particles in the system. This 
is the first time, up to the authors knowledge, that non-ergodic behavior in a system with 
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FIG. 4. Moments of the second moment (t 2 ) for the equilibrium and waterbag distributions. 



f(t) 




FIG. 5. Log-Log plot of the tail of the distribution function of sojourn times t s for the equilibrium 
case for TV = 50,000. The straight is a power law with a = 0.5849. 

long-range interactions is shown to be related to Levy nights in momentum space. We also 
performed a similar study for non-homogeneous states with similar results. A more thorough 
analysis for the HMF model and other models with long-range interactions is under course 
and will be the subject of a future publication. 
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